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A review is given on some recent developments in the theory of the Ising model in a random field. 
This model is a good representation of a large number of impure materials. After a short repetition 
of earlier arguments, which prove the absence of ferromagnetic order in d < 2 space dimensions for 
uncorrelated random fields, we consider different random field correlations and in particular the 
generation of uncorrelated from anti-correlated random fields by thermal fluctuations. In discussing 
the phase transition, we consider the transition to be characterized by a divergent correlation 
length and compare the critical exponents obtained from various methods (real space RNG, Monte 
Carlo calculations, weighted mean field theory etc.). The ferromagnetic transition is believed to 
be preceded by a spin glass transition which manifests itself by replica symmetry breaking. In 
the discussion of dynamical properties, we concentrate mainly on the zero temperature depinning 
transition of a domain wall, which represents a critical point far from equilibrium with new scaling 
relations and critical exponents. 



1 Introduction 

The Lenz-Ising model is probably the oldest and most simple non-trivial model for coop- 
erative behavior which shows spontaneous symmetry breaking EJI3. It has a vast number of 
applications ranging from solid state physics to biology. Its Hamiltonian is given by 

H = —'^JijSiSj — ^^HiSi (1) 

i,j i 

where Si = ±1. In a more restricted sense, the Ising model is understood to have cou- 
pling constants which are translationally invariant, J^ = J(|rj — rj\), where rj, rj denote 
points of a regular <i-dimensional lattice. The external field Hi is typically homogeneous, 
Hi = H{ri) = H\/i, or depends only smoothly on rj. Even under this restricted conditions 
the Ising model exhibits a remarkable complexity. In particular, if competing interactions 
are taken into account, the Ising model may show a large number of commensurate and 
incommensurate phases of modulated magnetization H. In reality, systems are never com- 
pletely translationally invariant. Compositional disorder, impurities and vacancies, lattice 
dislocations etc. lead to modifications in the Hamiltonian, which in many cases may be 
characterized by changes J(|rj — r^- 1 ) — > J(|rj — rj\) + 5Jij and H — > H + hi in TL. SJij 
and hi are no longer translationally invariant, but random quantities, characterized by their 
probability distributions. Typically the values of 5Jij and hi have a zero average and are 
uncorrelated for different bonds and sites, respectively. Let us briefly consider the different 
limiting cases separately. 

For hi = and 5 Jij <C Jij, we expect, that the influence of the randomness is not very 
dramatic. Higher order commensurate phases may disappear §. In the case of second order 
transitions the critical exponents can be changed H. First order transitions on the other 
hand may become second order Q. 

If SJij ^> Jij however, ferromagnetic order will in general be destroyed. A low tem- 
perature spin glass phase with < Si > = but < Si > 2 = qEA > may occur 0. Here 
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< . . . > denotes the thermodynamic average and the overbar the average over the disorder 
configurations. The latter replaces the spatial average over the sample. 

In this article we will consider the complementary case 5Jij = and hi 7^ 0, i.e. the 
Ising model in a random field. If not stated otherwise we will assume below 

hi = , h{hj = h 2 Sij and h -C J. (2) 

The exchange constant Jjj = J is assumed to be short ranged and non-zero between nearest 
neigbors < i,j > only. An alternative soft -spin description of the random field Ising model 
is given by the Ginzburg-Landau-Hamiltonian 

n GL = J d d r S.~a<f> 2 + i(V0) 2 + )u<f - h(r)d>(r)\ (3) 

with -oo < 0(r) < oo, a = a'{T - T c (0)), T c (h = 0) = 0(J),u > and 

h(r) = 0, h(r)h(r>) = h 2 5(r-r'). (4) 

Since its seminal discussion by Imry and Mai in 1975, this model is under intensive 
investigation both experimentally and theoretically. Results obtained before 1991 are sum- 
marized to a some extend in earlier reviews i>00lll. In the present paper we will therefore 
mainly concentrate on more recent findings and refer to the earlier work only whenever 
needed to make the paper more self-contained. 

As has been discussed earlier, the random field Ising model has a number of interesting 
realizations in nature. The most studied experimental systems are diluted antiferromagnets 
in a homogeneous external field where the combination of dilution and external field leads 
to a random field like effect for the staggered magnetization EL3. 

Another example is adsorbed mono-layers on impure substrates, e.g. Xe on a Cu(llO) 
surface. Here the mono-layer has two (or more) ground-states on the substrate. If one of the 
substrate lattice sites is occupied by a frozen-in impurity, it prevents additional occupation 
of this site with an ad-atom. Thus it acts locally as a symmetry breaking fieldlij. 

Further realizations of random field systems are binary liquids in porous media0, mixed 
Jahn-Teller systems EJ, diluted frustrated antiferromagnets EJ, hydrogen in metals and 
mixed crystals undergoing structural or ferroelectric transitionslli. Recently, an application 
of the random field Ising model for the understanding of cooperativity of protein folding 
has been discussed by Shakhnovich and coworkers B9 . Another more recent development 
is the discussion of the Anderson-Mott transition of disordered interacting electrons as a 
random field problemcil. 

The rest of the paper is organized as follows: In Section 2 we consider the stability 
of the ferromagnetically ordered phase, in particular in the presence of field correlations 
which deviate from Eq. ||. Section 3 is devoted to the discussion of the critical behavior. In 
Section 4 we discuss some dynamical properties of the model. Finally, Section 5 is reserved 
for miscellaneous subjects related to the random field problem. 
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2 Ordered Phases 



2.1 The stability of the ferromagnetic phase 

The pure Ising model Ho = — J J2<ij> Si^j is known to have a ferromagnetically ordered 
phase in all dimensions d > 1. An additional random field term 7i\ = —J2ihiSi will act 
against this order. When the random field strength h is sufficiently large co r n p ared to J, the 
system is disordered at low temperatures, as shown in an exact treatmentr 2 lr 3 l. The opposite 
case h <C J is of more interest. As it was first shown by Imry and Mala, the ferromagnetic 
ground-state becomes indeed unstable with respect to the formation of ill-oriented domains 
in all dimensions d < 2. 




+ 



Figure 1: domain of reversed spins 

If one reverses the spins (Fig. 1) within a domain of size R, the energy cost S ex is 
proportional to the domain wall area E ex ~ JR d_1 . For simplicity, we set the lattice 
constant equal to unity. This energy increase has to be compared with the energy gain from 
the interaction with the random field. Clearly, the average Zeeman energy vanishes for 
the ferromagnetic state. However, according to the central limit theorem, the mean square 
random field energy E^ F inside a region of volume R d is of the order h 2 R d . The energy Erf 
for a given region may be positive or negative with equal probability. It is always possible 
however to find a region enclosing an arbitrary point i such that -Erf > 0. Reversing the 
spins in this region yields an energy gain of —2Enp. The total energy of the domain of size 
R is therefore 

E(R) « JR^ 1 - hR d ' 2 . (5) 

For h <C J, E(R) is positive for d > 2, but negative for d < 2 if R is large enough. Thus 
the ferromagnetic state becomes unstable with respect to domain formation for d < 2. 
Later, Binder IHll was able to show, that in d = 2 dimensions the roughness of the domain 
surface generated by the random fields leads to an instability of the surface tension for 
R ~ R c ~ exp(J//i) 2 and hence there is no long range order also in d = 2 dimensions. 
Aizenman and Wehr@ proved indeed rigorouly, that in d < 2 dimension the random field 
produces a unique Gibbs state, i.e. the absence of any phase transition, in agreement 
with the naive expectation from the Imry-Ma argument. In higher dimensions the surface 
roughness has no influence on the long range order. 

So far we have assumed, that domains are compact and do not include smaller domains 
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of reversed spins. Moreover, entropic effects were neglected. Indeed, it could be shown, 
that 

(i) domains within domains, which would renormalize h — > h(R), are rare if d > 2 and 
can hence be neglected @, 

(ii) although there is a large number N(R) ~ exp (cR^ 1 ) of contours for a domain of 
size R, these enclose essentially the same volume and the same random fields. Here c is 
an unknown numerical coefficient. Thus, with probability one there is no contour which 
encloses a random field gain which is larger than the surface energy loss if d > 2c3, 

(iii) ,— thermal fluctuations are irrelevant at low temperatures and can hence be ne- 
glected c3. 

Thus, one concludes, that the lower critical dimensions for the ferromagnetically ordered 
phase is di = 2. 



2.2 Different field correlations 

So far we have assumed, that both the interaction between spins and the correlations be- 
tween the random fields are short range. In this paragraph we will briefly consider the 
extension of the Imry-Ma argument to long range interaction and correlations, respectively. 

Let us first consider the case of infinite-range interaction, Tto = —(J/N) J2i<j SiSj- 
Here N denotes the total number of spins in the system. The energy increase by flipping 
a domain is now of the order JR d and hence always larger than the energy gain from the 
random field. The system is therefore ordered at low temperatures in all dimensions. 

Another type of long range interaction comes from dipolar forces. For nearly spherical 
domains these lead to an additional term gR 2d ~ 3 in Eq. ||, where we assumed three dimen- 
sionaldipolar interaction (vanishing as i? -3 ) between dipoles arranged in a (f-dimensional 
spaceca. g is a measure of the strength of the dipolar interaction. It is easy to show, that 
this term does not change the lower critical dimension di = 2. In general, domains may not 
be spherical and indeed will tend to take a cigar-like shape with the long axis parallel to the 
spin direction in order to lower demagnetisation factors, but this will not change di either. 
A slightly different treatment is necessary for diluted antiferromagnets in a homogeneous 
external field, which are a good experimental realization of the random field Ising model, 
but the conclusion d\ = 2 applies also to this casecS. 

The situation becomes different if we allow long-range correlations of the random fields, 

i.e. 

H~h~= h 2 \ri -rj\ d '- d for i^ j and 7tf = h 2 . (6) 

For d! > the lower critical dimension changes to di = 2 + d'S. Examples are random fields 
correlated along lines (d' = 1) or planes (d' = 2). A special example is the d-dimensional 
quantum ferroelectric with uncorrelated random fields at T = 0, which can be described by 
a (d+ 1) = D-dimensional classical model and d! = 1, since the disorder is correlated in the 
(imaginary) time direction. Thus Dj = 2 + 1 = 3 and hence d\ = D\ — 1 = 2, i.e. quantum 
effects do not change di in this caseEll. 
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An opposite situation exists in systems with anti- correlated fields® 

w?=( -r i |-" ,+2 '" . \*i (7) 



where c{ji) is a constant determined by the anti-correlation condition J2j hihj = 0. If this 
condition were not satisfied, the properties would be the same as for uncorrelated fields. 
c(yu) vanishes as \x approaches zero. Therefore, the limit [i = corresponds to uncorrelated 
fields. 

Let us start with the case T = 0. In order to obtain d[ we have again to calculate 

£rf~ E W ~ h2 ( Rd ~ %11 + Rd ~ l ) (8) 

|ri- ri |<il 

where the second term on the rhs is a surface contribution. Comparing .Erf with the surface 
energy Ji? rf_1 we obtain d\ = 2(1 — /i) if \i < 1/2. For fi > 1/2 the random field does not 
affect the lower critical dimension and hence di = 1. 

However, as we will show now, these are not the lower critical dimensions at finite 
temperatures. To illustrate this point we consider random field dimers, which consist of 
anti-parallel random fields ±h on two neighboring lattice sites and which are a particular 
realization of anti-correlated random fields corresponding to large values of fi. Clearly, 
at T = 0, dimers (compare Fig. 2) do not change the energy E(M = ±1) of the two 
ferromagnetic ground states with magnetization M = ±1 and hence E(M = ±1) = Eq. 
Random field dimers act here similar to (oriented) random bonds, which can lower the 
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Figure 2: Different configurations of random fieid dimers 

surface energy of a domain, as follows from the second term in Eq. ||. 

The situation becomes however different if we consider the free energy difference of the 
states M ~ ±1 at non-zero temperatures < T <C h <C J- For a single isolated dimer 
(Fig 2. a) we find from the low temperature expansion of the free energy again no symmetry 
breaking effect 

F(M « ±1) = Eq- Te- 2zJ ' T (N-2 + e 2h ' T + e~ 2h l T ) + O ( e - 4zJ / T ) . (9) 

Here z denotes the number of nearest neighbors of a spin. The second term on the rhs of 
Eq. |j| arises from flipping a single spin only. On the other hand, for two dimers forming a 
pair as in Fig. 2b, c we get instead 

F(M « ±1) = E - Te~ 2zJ ' T (N - 3 + e^ h ' T + 2e ±2h ' T ) + O (e~ izJ / T ) . (10) 
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Thus 

F(l) - F(-l) « 2Te~ 2zJ/T (sinh (4/i/T) - 2 sinh (2/i/T)) = 2/i cfr (11) 

where /i e g can be considered as the effective random field strength of the dimer pair. Similar 
symmetry breaking contributions follow in higher order in exp (—2zJ/T) also from spatially 
separated dimers as in Fig. 2d. and from other realizations of anti-correlated fieldstH Thus 
anti-correlated random fields generate uncorrelated random fields, as can alternatively be 
seen from a perturbative treatment of the Ginzburg-Landau-Hamiltonian@. 

In d = 2 dimensions we come therefore to the surprising conclusion, that although a 
system with anti-correlated fields has a ferromagnetically ordered ground state, ferromag- 
netic order is destroyed at all non-zero temperatures. The correlation length diverges for 
T — > as exp ( J/h e g) 2 where h e s is given in Eq. 11. We note finally, that anti-correlated 



random fields may play a role in the description of binary fluids in VycorS. 

2.3 Interface properties 

In the argument of Imry and Mai one assumes, that domain walls are smooth such that 
their surface area is proportional to R d_1 if R denotes linear extension of the domain. A 
closer inspection shows, that this is indeed true in d > 2 dimensions. For simplicity we 
consider here a domain wall with a vanishing mean curvature, which can be considered 
to be a small part of the surface of a large domain. If we denote the interface profile by 
z(x) where x denotes a (d — l)-dimensional vector parallel to the mean interface plane and 
r = (x_, z) , the interface Hamiltonian is 

Hi = Jd^x |a + ^r(Vz) 2 -2M^ 2( " ) dz / /i(x,z / )|- (12) 
Here ctq and T denote the bare surface tension and the wall stiffness, respectively. 



Eq. 12 can be derived in principle from the lattice Hamiltonian of the random field 
Ising model. <jq and V depend then in general on J, h and T in a complicated manner. 
At T = 0, cr ~ J and V ~ J(J/h) 2 ^ 3 ~ d ^ for d < 3 and T ~ J e c ^ J / h ^ for d = 3 H. 
For 3 < d < 5 dimensions and h < h c domain walls are flat because of the influence of the 
periodic potential of the underlying lattice, and hence cannot be described by Eq. IT! 



the wall undergoes a roughening transition such that for h > h c Eq. |l2j again appliesl 33 H 34 » 
In the following we will assume, that we are always at h > h c . 

The properties of the random field h(r) are still given by Eq. [|. The smoothness of the 
domain wall can be checked by considering the interface roughness 

w(L) = <{z{x l )-z{x 2 )Y > 1/2 ~ L<. (13) 

\x.i-x 2 \=L 

If the exponent £ is nonzero, but less than one, the interface is termed self-affine. 

A simple Imry-Ma argument for the energy of a hump of a linear extension L and height 
w gives 

^f^yvw^- (14) 
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The expression on the rhs of Eq. 14 defines the Larkin length L c , 
This result Eq. 14 has been confirmed in d 

and in d 



Section 4.1 



to which we come back in 
5 — e dimensions by a functional 
2 by a mapping on the Burger's equationEEl. 



renormalization group calculation 
The roughness exponent Q = (5 — d)/3 is smaller than one in all dimensions 2 < d < 5. For 
d > 5 dimensions domain walls are always flat (one could also say, that h c — > oo). Thus 
in all dimensions d > 2, the typical gradients of the walls are small ~ L^ 1 ~ L( 2 " d )/ 3 and 
hence walls are smooth on large scales. On the other hand, in d = 2 dimensions domain 
walls are rough on all scale, which leads to a vanishing total surface tension on length scales 
L > R c ~ exp(ca r 1 /3/ l -4/3) 

~ exp [c'(J/h) 2 } c3. c and d are constants of order one. 



3 Phase Transition 

3. 1 Order of the transition 

In the last section we convinced ourselves, that the Ising model in a short range correlated 
random field has an ordered phase in more than two space dimensions. A natural question 
is then about the order of the transition. 

The mean field approximation, which is exact for the infinite-range model, predicts a 
second order transition for a Gaussian distribution of the random field strength. For a 
bimodal distribution the transition however becomes first order if h is larger than a critical 
valued. For a model with short range interaction earlier numerical simulations^ and high 
temperature expansionsed seem to confirm a first order transition, but now independent of 
the field distribution. 

More recent Monte Carlo simulations for bimodal and Gaussian random field distribu- 
tions, however, do not show a latent heat EH or phase coexistence E3, but a jump in the 
magnetization M. The latter may be related to a logarithmic dependence of M on the 
reduced temperature. Also recent high temperature series expansion up to 15 term s sh ows 
a continuous transition for both distributions e3. On the other hand, Swift et al. find 
in a recent study from an exact determination of the ground states in higher dimensions 
{d = 4) a discontinuous transition for a bimodal field distribution whereas for a Gaussian 
distribution and in d = 3 dimensions the transition is continuous. 

Although there is so far no proof for a continuous transition in neither case, we will 
assume in the following that the transition is indeed continuous, which seems to be at least 
true for a Gaussian random field distribution. 



3.2 Scaling laws 

Under which conditions are random fields relevant at a critical point? To answer this ques- 
tion we consider the system at the true transition temperature T c {h) in a finite homogeneous 
field H. We divide the system into blocks of linear size £. Here £ denotes the correlation 
length which is related to H by H ~ J£~ /35 /* y . Here we follow the standard notation for 
the critical exponents a, j3, 7, 6, v etc. ea. The random field produces an additional excess 
field 5H ~ cht^ d l 2 where c is a constant of order one and arbitrary sign. Approaching the 
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critical point H — > we expect a sharp transition, if for £ — ► oo, 5H(£)/ H(£) — > 0, i.e. for 

d #5 2 + 7-g 

2 " V = —2^' (15) 

On the rhs of Eq. ^ we used a (i-independent scaling law which is assumed to be valid also 
for random field models. With the mean-field exponents (3 = v = 1/2, S = 3 we find, that 
random field effects are negligible only for d > d c = 6. Thus random fields are relevant for 
d < 6 cB. For d < 4, where the pure model has already non-classical exponents 
etc. we get from eq. Eq. |l^, and using the conventional hyper-scaling law v^'d = 2 — a(°) 
the relation 7^) < 0. Since the susceptibility exponent 7^) has to be positive, we have 
to conclude, that the random field is also a relevant perturbation if we use the critical 
exponents of the pure model and moreover, that the conventional hyperscaling breaks down 
in random field systems. 

The condition for the applicability of the linear response theory, 5H(£) < H(£), is 
violated if 

/ t\ 2v/(2/3S-ud) 

t>tu>~ J • P*> 



Eq. 16 represents the Levanyuk-Ginzburg critical region, where random field effects cannot 
be ignored, since different volumes perceive different effective fields H + 5H. The resulting 
transition could be 

(i) smeared (as one could naively expect) or 

(ii) 1st order (as ruled out in the previous paragraph), or 

(hi) 2nd order with modified exponents and modified scaling relations (because we would get 
the unphysical result 7 < with conventional hyper-scaling relation). 

Let iL&_assume that case (iii) (as simulations suggest) and the modified hyper-scaling 



relation 
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u(d-6) = 2-a (17) 
with some new exponent 9 apply. Then relation Eq. ^ leads to the condition 

#>7/V = 2-7] (18) 

an inequality which has been proven exactly by Schwartz and Soffer for a number of random 
field systems 111. 



A heuristic argument to make Eq. |17| obvious follows the original idea of Pippard: 
the free energy of a correlated region Ft scales as the energy, which is necessary to flip 
the region. In pure systems this energy is ksT whereas in random field systems it can be 
estimated as h£ e . Our previous Imry-Ma argument would give 6 ~ d/2, but we have to 
take into account, that the random field gets renormalized close to criticality. Thus 

F f w j^-( 2 - Q )/^ „ k B T + hi e (19) 
from which we can immediately read off Eq. [l7]. is a third independent critical exponent 



that describes the irrelevance of temperature as is plausible already from Eq. 19 
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The appearance of a third exponent is related to the existence of two different correlation 
functions which scale independently. In particular, at T = T c (h) 



and 



G(k) = (S k S_ k ) - (5 k ) (5_ k ) ~ AT 2+7 > 



C dis (k) = (S k ) (5_ k 
where »§ k = h= X^exp (ikr^)Sj andBB 



1] = 2 + T] 



(20) 
(21) 

(22) 



Note, that Cdi S (k) vanishes for pure systems. 

Besides the Schwartz-Soffer inequality there is a number of further inequalities. From 
(3 > follows with Eq. 17 and Eq. 22 f\ > 4 — d. Further inequalities are 7 > E§ and 



d—X > 9 cH. The latter inequality follows from the fact, that the domain wall energy scales 
as ^~^l u with \i = v(d — 1 — 6) > 0. Thus, the values of j3/v and 9 are restricted to a 



pentagon (see Fig. 3). Schwartz and Soffer have further claimed that Eq. 18 is fulfilled as 
an equality^!. Finally, there is the Harris- inequality v > 2/d. 



d=2+e 




9/d 



Figure 3: The domain of allowed values of ijh = d — /3/v and y.j = 0. The arrow denotes the result of the 
perturbative expansion in d — 6 — e dimensions. Most of the results for the exponents found numerically are 

located in the upper left corner 



3.3 Renormalization group 



Since 9 > 0, it is plausible from Eq. 19, that in the critical region thermal fluctuations 
are less relevant than those of the frozen-in disorder. In a renormalization-group (RNG) 
treatment this feature is reflected in the existence of a T = fixed point, which is believed 
to describe the critical behavior of the random field Ising model. 

Although there is so far no satisfying RNG analysis, we present here briefly a rough 
sketch of it, assuming a continuous transition up to T = 0. We follow thereby closely the 
earlier work of Bray and Moore Ej. We start with the observation, that the free energy 
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density / can be written in the form / = Jf(T/J,h/J,H/J). Let us imagine to carry 
out the RNG coarse-graining transformation, with length scale factor b, corresponding to a 
reduction in the number of degrees of freedom by a factor b d . The transformation generates a 
flow in the space of the naive scaling fields T/J, h/ J, and H/ J, which eventually terminates 
in one of the fixed points of the system. The existence of three fixed point will be assumed 
(Fig. 4), in addition to the trivial, high temperature fixed point: 

(i) A totally unstable "thermal" fixed point C at T = T c , H = h = (the random field is a 
relevant perturbation, see our discussion in |3.2| ) . 

(ii) A fixed point R at T = H = and h = Iir which is unstable in two, but stable in one 
directions and is therefore a critical point. 

(hi) A totally stable fixed point FatT = h = H = 0, which corresponds to the low 
temperature phase for d > d\. 




Figure 4: Schematic renormalization group flow of the random field fsing model 



In general the RNG procedure generates also new terms in the Hamiltonian. We will 
assume, that these terms are irrelevant in the RNG-sense and can therefore be neglected. 

In order to calculate the critical behavior we have to linearize the RNG flow close to 
the fixed point R. The eigenvalues and eigenvectors of the linearized RNG-transformation 
deliver the critical exponents and scaling fields. Phenomenological arguments concerning 
the RNG flow suggest 

j, t = j(h - h R ) + cj and - (23) 

as the scaling fields. Close to the fixed point R J, r and H transform under the RNG 
coarse graining as 

J -» J' = Jb yj 

t^t' = rby- (24) 
H -» H' = Hb VH . 
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From Eq. 23, Eq. 24 and the invariance of the partition function under the RNG 
transformation we get for the singular part of the free energy density f = J • f 



/( J' T 'J 



bVJ- d f f-b- yj , Tb VT , —b VH - yj 

V J J 



and similar for the correlation length £ 



4 (j> T >j)= b z(j b ~ VJ > TbVT >J bVH ~ yj 



(25) 



(26) 



The critical exponents follow from Eq. ^ and Eq. |26| in the usual way, their relations to 
eigenvalue exponents yi are summarized in the following list: 



j3 = v(d- y H ) v 

7 = v(2y H -yj-d) rj 

8 = (vh - y,j)/{d- y H ) fj 



2 + d + yj - 2y H 
A + d-2y H 



= yj 



(27) 



These scaling relations are in agreement with Eq. [l?] and Eq. 22. The existence of a 
zero-temperature fixed point corresponds to a positive value of yj = 6. Alternatively, 
— yj can be interpreted as the eigenvalue exponent of the temperature or, by simple scale 
transformation <\> = y/Tc/)' in Eq. |3[ of the coupling constant u. With this interpretation 
of —yj, the modified hyper-scaling law Eq. 17 has been first derived by Grinstein §. The 
RNG program scetched above has been performed in real space using Migdal-Kadanoff or 
related approximations by a number of groups EHilijl'H (see Table IT]). 



3.4 Dimensional reduction and weighted mean field approximation 

A convenient starting point to study the critical behavior is the Ginzburg-Landau Hamil- 
tonian H,QL{cf)}. Including a random field term, Young was able to show that the most 
singular terms in the perturbation theory for this model follow then from tree diagrams, 
which can alternatively be obtained from an iterative solution of the saddle point equation 



5H G L/6<j> = -V4> 2 + V\4>) - h(r) 



0. 



(28) 



Here V((f>) = ^dcf) + ju<p . The contributions from tree diagrams lead to an exact relation 

between the critical exponents (d) of the random field system in d-dimension and 
those of the pure model in {d — 2)-dimensions 



\f F \d)=Xr(d-2) 



(0), 



(29) 



the so called dimensional reduction. This equality has been proven, starting from Eq. |28| , 
also in a non-perturbative wayE3. Comparing Eq. with eq. Eq. 17 yields 8 = 2 in all 



dimensions. Since the lower critical dimension of the pure Ising model is di = 1 one would 
conclude d) = 3, in disagreement with our findings from the previous Section. 

Apparently, perturbation theory is inappropriate to deal with this type of disorder which 
is characterized by a large number of local minima in the energy landscape. Indeed, since 
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the perturbation theory can be generated from Eq. 28, all saddle point solutions enter 
expectation values of physical quantities with the weight ±1, which is clearly the wrong 
way. 

More recently Lancaster et al. H proposed a weighted mean field theory which takes 
into account all solutions mf of the mean field equations 

dF/drrij = — ,' ^rrij + Tarctanhm^ — hj = 0, i = l,...,N. (30) 

Here Ej denotes the sum over the nearest neighbors to the site i. hi is taken from 
a bimodal distribution. Thermodynamic quantities are then calculated as a sum over all 
mean field solutions a with the Boltzmann weight exp (— F a /T)/ Ea ex P ( — F a /T) where 
F a = E a — TS a is the mean field free energy and 

E a = -\J E mfmf-Ehirnf 

<ij> i 

• (31) 

Sa = ~ E i {(1 + mf) ln (! + m i) + (! - m i) ln (1 - m i) - 2 ln2} 



We note, that Eq. pQ is a lattice version of Eq. 28 with V(rrii) = —T c m,{ + Tarctanhm. 



and T c = 2d J. Since Eq. |30| represents the saddle points of F a and thermal fluctuations are 
expected to be irrelevant for the critical behavior, in principle one should be able to calculate 
the true critical exponents from this approach (at least at sufficiently low temperatures). 
Starting from high temperatures T > 7^ pure ) there is typically only one solution to the mean 
field equation. Decreasing T below a temperature T* > T c , the number of solutions starts 
to grow rapidly. For a 32 3 periodic lattice and an accuracy \mf — mj\ ~ 1CP 3 , this number 
becomes of the order 250, although the number of those with a large weight increases more 
slowly. The critical exponents obtained in this way are summarized in Table |I|. Lancaster 
et al.Ea identify T* with the temperature Trsb, at which replica symmetry breaking occurs 
(see Section 3.6). 

3.5 Two or three independent exponents ? 

In contrast to conventional critical points, which are characterized by two independent crit- 
ical exponents, our schematic renormalization group calculation suggests, that the random 
field Ising model is characterized by three independent exponents. 

In an early publication, Aharony, Imry and MaS suggested the existence of an exact 
exponent relation 

(9 = 2-7? (32) 

which implies 2ry = f\ and reduces the number of independent exponents again to two. 
Relation Eq. ^ can indeed be made plausible by estimating the free energy of a correlated 
droplet close to the critical point. With M(r) « xh^ d ^ 2 ~ h^ ri ^ d l 2 for the local 



magnetization, where x = G(0)/T denotes the susceptibility, and Eq. |19| we get 

d d r fc(r)M(r) - h 2 ^ ~ f (33) 



d 



12 



Table 1: Critical exponents. 



Reference 


= yj 


v = l/y T 


(3 = v{d-y H ) 


2rj-fj = 6-2 + r] 


MIA 


a — 4 


1/2 


1/2 


a — 6 


PT b d = 6-e 


2 


i/(°)(d-2) 


/3(°)(rf-2) 


?/ )(d-2) 


d-> 2 


1 





2 





weighted MFA@ 


1.51 


1.25 ±0.11 






MC c (d=3) 

Rieger and Young Cj 
RiegerFH 


1.56 ±0.1 
1.53 ±0.1 


1.6 ±0.3 
1.1 ±0.2 


0.003 ±0.05 
0.00 ±0.05 


0.12 ±0.12 
-0.03 ±0.15 


HTSE0 

Gofman et al. c3 


(7 = 2.1 ±0.2) 









RealspaceJING 

Dayan et al. E3 . . 

Newmann et al. El 
Cao and Machi-aL_l 
Falicov et al.Ej 


1.56 

1.00 ±0.05 
1.5 

1.4916 ±0.0003 


1.39 

1.49 ±0.008 

v = 2.25 

i/ = 2.25 ±0.01 


W 

1.66 ±0.01 
0.02 

0.0200 ±0.0005 


0.12 
1.24 
0.02 

0.001 ±0.001 



which gives Eq. 32. Later Schwartz et al. cEI have claimed, that there is an exact proof for 
this relation. However, all these approaches use in one or the other way linear response 
arguments and have therefore to be considered with caution. Numerical calculations show, 
that Eq. ^ is indeed fulfilled within the accuracy of the calculation. 

Since the numerical determination of exponents in random systems is typically hampered 



by the existence of considerable error bars, which makes a confirmation of Eq. 32 difficult, 
Gofman et al. El considered the even stronger relation 



A 



lim 



T 2 C dls (0) 



(34) 



h 2 G 2 (0) 

which should be fulfilled according to El. 

The exponent scaling gives for Cdj s (0)/G 2 (0) ~ £ 2r )-v^ which would diverge (compare 
Eq. 18 and Eq. [2?]) unless Eq. ^ is valid. Gofman et al. determine A from a 15 terms 
high temperature series expansion for G(0) + Crfj s (0) and G(0) in d = 3, 4, 5 dimensions 
for different values of the random field strength h. In particular, for d = 3, A = 1 ± 0.003 

(but not a proof!). Thus, despite 



which is an impressive confirmation of Eq. 34 and Eq. |3 
of all efforts to prove Eq. this problem has still to be considered as unsolved. 



a MFA: mean field approximation 

6 PT: perturbation theory 

C MC: Monte Carlo simulation 

d HTSE: high temperature series expansion 
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3.6 Replica symmetry breaking (RSB) 

The disorder average can conveniently be performed using the replica trick. E.g. the average 
free energy can be written as 

F = -TWZ = -T lim - fZ^-1) (35) 

n->0 n 



where Z n = (Tr exp (-H/T)) n = Tr exp (-(l/T)W n {5f }) is the conventional partition 
function of the translationally invariant replica Hamiltonian 

H n {Sf } = - E E ( J v 6 <*P + If • ( 36 ) 

a,/3 i,j \ / 

An alternative formulation uses the Ginzburg-Landau model as the starting point with the 
corresponding replica Hamiltonian 

n GL , n {r(r)} = E H GL m - E / d d r^-r^- (37) 

Mezard and Young use an m-component generalization of Eq. [j?], <p a (r) — ► (j) a {v) = 
(4>f(r),i = 1 . . . m), to determine the correlation functions G a @(k) and C^f(k) from Bray's 
self-consistent screening approximation (SCSA). Note, that this generalization changes the 
lower critical dimension from di = 2 to di = 4, which follows both from perturbation theory 
and, because of the continuous symmetry of the order parameter, also from the Imry-Ma 
arguments. However, dimensional reduction is expected to break down also in this case £3. 

The SCSA is a truncation of Dyson's equation which is exact to order 1/m. Assuming 
a replica-symmetric solution, Mezard and Young find the dimensional reduction result 
77(d) = fj(d) = r]( \d — 2) = 0(l/m). This replica symmetric solution is however unstable 
with respect to RSB. Using the hierarchical RSB scheme of Parisi, they find 77 = rj(°\d — 2) 
to order 1/m, but the value fj is now altered. The difference 277 — 77 = c/m is small 
for any m. The constant c can in principle be determined from the set of self-consistent 
equations. This calculation has been extended by Mezard and MonassonS and de Dominicis 
et al.EEl, who determined the temperature Trsb > XL where the replica symmetric solution 
becomes unstable. In particular, de Dominicis et al.tJ consider a more general coupling term 
^ k A^(k)^(^ k between the replicas in Eq. 37 and determine Trsb from the divergence 



of the A-susceptibility d 2 Z n / dA a p(p)dA^s(p') (to reach this goal, in practice they use the 
Legendre-transform of Z n ), which is related to the standard spin-glass susceptibility. It 
turns out, that the spin-glass transition (in the sense of an Almeida-Thouless line), which 
is believed to take place at ?rsb, always precedes the ferromagnetic transition (see Fig. 5). 
For d close to six dimensions they get in particular 

\TBSB-T c (h)\ « J{h/Jfl^ d ). (38) 



This has to be compared with the Levanyuk-Ginzburg criterion, Eq.16, which gives a ran- 
dom field controlled critical region of size 

|T LG - T c (h)\ « J(h/Jfl^ » ITrsb - T c (h)\. (39) 
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For T c < T < Tlg we expect, that the saddle point equations Eq. Eq. 30, will have 
several solutions, which signals the failure of perturbation and linear response theory. It 
is at present unclear, whether RSB (and hence the breakdown of dimensional reduction) 
occurs in the whole random field critical region, as one would naively expect and as was 
found in a more recent study by Dotsenko and MezardEa, or, as Eq. [3^ suggests, is restricted 
to the much smaller temperature interval around T c (h) (h -C J). We note finally, that the 
spin-glass order parameter q^A 1S non-zero at all temperatures. In particular, outside the 
critical region where linear response theory applies < Si > 2 = qEA ~ J2ij XiiXijhihj = 
h 2 J2iXh > 0- Here T%ij =< SiSj > — < Si >< Sj > denotes the Fourier transform of 
G(k). 

Although we have here only discussed the paramagnetic phase, one expects a similar 
behavior if one reaches the ferromagnetic transition line from below. 



4 Dynamical Properties 

4-1 Zero temperature interface depinning 

At finite temperatures the two phases with up and down magnetization can coexist only 
for vanishing strength of the uniform external field H = 0. This is not the case at zero 
temperature, where the disorder (and possibly also lattice effects) lead to a pinning of the 
wall separating the two domains. In order to get the wall depinned, the external field has 
to overcome a threshold H c . Thus, the coexistence surface consists of two parts, for < 
T < T c {h) it is restricted to H = 0, whereas for T = it is given by —H c (h) < H < H c (h) 
(see Fig. 5). 




Figure 5: Coexistence surfaces of the random field Ising model 

In this section we consider the behavior of the interface in the vicinity of the depinning 
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threshold and determine H c . The equation of motion of an over-damped wall isLJ 

A— = -dHr/dz = TV 2 z + 2M{H + h(x, z)). (40) 

This equation is highly non-linear because of the last term on the rhs. A denotes the inverse 
mobility of the interface. We assume here, that the interface motion is over-damped. We 
will show later, that an inertial term is indeed irrelevant close to the depinning threshold. 
It turns however out, that it is important to assume a finite correlation length a for the 
random field in the ^-direction (which is the direction of motion of the interface), i.e. we 
replace Eq. || by 

h(x, z)h{x', z<) = h 2 5 d -\x - x!)A(z - z') (41) 

with /_ A(z) = 1. a is of the order of the lattice spacing, i.e. of order one in our units. 
A(z) = — A(— z) is a monotonically decreasing function of z for z > and decays rapidly 
to zero over a finite distance a. The width of the random field correlator perpendicular to 
z turns out to be irrelevant, we assume therefore, that it is smaller than any other length 
of the problem. 

The physics of the interface close to the depinning transition is characterized by two 
emerging important length scales, the Larkin length L c and the dynamical correlation length 
£„. For weak disorder, J 3> h, L c is the length scale on which the typical distortion of the 
interface is of the order a. We have to compare here the curvature force FL d ~ 3 a with the 
random force ~ M/iA 1 / 2 (a)L (d_1)/2 . In the case of weak disorder, which we consider here, 
Ta 3> M/iA 1 / 2 , and the random force wins over the curvature force only on length scales 

L > L c « [(Ta/Mhf/Aia)} 1 /^. (42) 

and d < 5. For L < L c (or d > 5), the interface is flat and hence cannot be pinned, since 
the total pinning force ~ L^" 1 )/ 2 is always smaller than the driving force ~ L( d_1 ). On the 
other hand, the interface is able to explore the inhomogeneous force field on larger length 
scales L > L c . It follows, that the maximum pinning force on a piece of interface of linear 
dimension L > L c is of the order (L/L c )^ d_1 ^M/i[A(a)Lc d 1 ^] 1/ ' 2 , which leads to a threshold 
field of the order 



M c M 



fhM\ 2 . 

br' A(oj 



2/(5-d) 

(43) 



For small h, H c ~ /i 4// ( 5 ~ d ), i.e. H c ~ h 2 in d = 3 dimensions. Close to the critical point 
h ~ h,R, H c ~ r^ _1 /M ~ t@ s . Note however, that T carries its own /i-dependence if Eq. 40 
is derived from the bulk Hamiltonian (compare Section 2.3). Then for h — ► 0, H c becomes 
exponentially small in d = 3 dimensions. 

To understand the nature of , it is convenient to describe the interface in a co-moving 
frame z(x, t) = h(x, t) + vt where v =< z > denotes the mean interface velocity. Hence 

dh 

A— = TV 2 h + 2MH -Xv + 2Mh(x, vt + h). (44) 
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Let us consider the motion of a domain wall over some typical obstacle formed by a random 
field cluster, which is assumed to be hit by the wall at t = 0. For small t, the wall will 
locally be stopped whereas other parts will continue to move forward by a distance vt. If the 
moving wall behaves self-afnne under the action of the random field (an assumption which 
we will prove later), it will form locally a bump of typical height h ~ t^ z , where £ and z 
are the non-equilibrium roughness and the dynamical index of the interface, respectively. 
Since C, < z, on time scales t > i„ ~ v ~ z /( z ~0 we fi nc [ vt ^> h and hence the non-linearity 



in Eq. 44 can be neglected. On the time scale t v the local perturbation by the obstacle 



spreads out on the interface over a length scale 

& ~ t\j~ z ~ u-V^-C) (45) 

which has apparently the meaning of a dynamical correlation length. Thus, on time and 
length scales t S> t t , and L > ^, the effective interface equation is the linear Edwards- 
Wilkinson-equation 

d z 

Aeff^ = r eff V 2 z + 2M{H eS + h eR (x, vt)) (46) 

but with renormalized parameters. In particular, we will find H c q = H — H c . In the rest of 
this section we determine the dependence of the effective parameters on (H — H c ) due to 
the non-linear effects on scales L <C £ v . Once these dependencies are known, Eq. ^ can 
be treated easily, it yields in particular in d = 3 dimensions a logarithmic roughness of the 
interface for L>(„. 

We start with the derivation of some scaling relations. If the transition from the moving 
to the pinned state of the interface is continuous, as we will show below, we expect a power 
law behavior of the velocity 

v~(H- H c f (47) 

and hence ~ (H — H c )~^^"~^' . One condition for a continuous transition follows from 
a Harris-like argument: the fluctuations 5H(£ V ) of the threshold force in the correlated 
volume ^ _1+< ' occupied by the rough interface should be small compared to (H — H c ) 0. 
Thus, (H - H c ) > 6H e (£ v ) ~ f-^+ft/ 2 ~ (H - tf c )*(<*-l+C)/2 or 

~ * >_—^ (48) 



z-C d-l+C 

where we introduced v via £ v ~ (H — H c )~ u . 

For large velocities the non-linear terms of Eq. |44| can be treated by standard pertur- 
bation theory l§. For the renormalized inverse mobility we obtain e. g. 

Aefr = A (l - Cl (£(°)/L c ) 5 - d + c 2 (^/L c ) 2 ^ - ...) . (49) 

Here L c is given by Eq. with A (a) replaced by |A"(0)|a 2 and = (Ta/vX) 1 / 2 is the bare 
dynamical correlation length. Similar expressions can be obtained for the renormalization 
of A"(0). There is no renormalization of T, however, because of an exact tilt symmetry 
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The perturbation theory clearly diverges for v — > 0. It turns out, that a functional 
renormalization group calculation is required to find the true critical behavior since A"(0) 
develops a singularity under the renormalization group transformation if b ~ L c . The 
resulting functional flow equation for the random-field correlator is 



dA(z) 
dlnfe 



(Z-2QA(z)+SzA'(z)-c^ 



z) - A(z)A(0) 



(50) 



where e = 5 — d. The fixed point function has a cusp-like singularity at the origin A*(z) = 
A(0) + a\\z\ + d2Z 2 + . . . where a\ and a2 are of the order e. It is important to note, 



that despite of the fact, that we were able to obtained H c in Eq. (43| from a simple scaling 
argument, the determination of H c from a straightforward calculation requires the existence 
of a non-zero value of A'(z — > + ) ~ H c . The roughness exponent £ follows from the 
fixed point condition to be £ = e/3, which is identical with the equilibrium roughness 
exponent. The correlator of the effective random field h e ff, appearing in Eq. 46 is A e ff(z) ~ 

(VCt,) (g - 2f) A*(z(L c /C„)f). 

The renormalization of the mobility is coupled to the renormalization of A e ff(z) and 
yields A e // ~ (L c /£„) 2e//9 . Together with the scaling relations S 



(51) 



which is valid to all orders in e E3, this gives for the exponents to order e = 5 — d 



c 



2 

2- -e, 
9 



(52) 



Narayan and Fisher have claimed, that ( is correct to all orders in e, but numerical 
calculations show deviations from this result in low dimensions^. With these replacements, 
all terms in the Edwards- Wilkinson-equation Eq. ^6| scale as H — H c . It is easy to see, that 
at the depinning fixed point an inertial term is irrelevant since it scales as (H — H c ) x with 



A = l + 2§Ej > 1. 



Thermal fluctuations will smear out the transition such that the velocity is non-zero 
for all driving forces HM. For H ~ H c 



v m T&p^{{H - H c ) /T l /P) (53) 

with ip(x) ~ for x — > oo and tp(x) ~ const for x — > 0. The value of /? is at present 
unknown. For H <C H c , the velocity is exponentially small and results from a creep motion 
over the barriers formed by the random field 

v ~ exp (—CMh 2 /HT) (54) 

with C = 0(1) 0. 
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4-2 Critical dynamics 

Villain @ and D.S. Fisher have proposed to consider the spin dynamics close to T c (h) as 
due to a kind of domain reversal, where now domains within domains should be taken into 
account through the exponent 6 introduced in Eq. [l7] and Eq. [[{]. The reversal of domain 
of size £ (the correlation length) is associated to an energy barrier of order h£ e , and the 
Arrhenius law gives the following expression of the relaxation time 



i re i exp( const. 



T-T r . 



V0 



)• (55) 



This formula is in contrast with usual power laws and reminiscent of the Vogel-Fulcher 
law observed in amorphous systems O. We note, that conventional perturbation theory 
gives t ve \ ~ \T — T c \~ uz with dynamical critical exponent z = 2 + 2n E3. Because of the 
exponential increase of the relaxation time, random field systems will rapidly fall out of 
equlibrium by approaching the transition. 

4-3 Metastable domains 



According to our finding of the Sec. 2.1, the 3-dimensional random field Ising model should 
exhibit ferromagnetic order at low temperatures. Experiments on random field systems 
are mainly perfor med with diluted antiferromagnets in external field and show pronounced 
hysteresis effects II30>EJ. In particular, no long-range order is found for 3d systems cooled 
from the high-temperature phase in non-zero field H. Neutron scattering experiments in 
three dimensions yield a non-equilibrium correlation length 

Rc(t) ~ H~ UH ■ r(t) with v H « 2.1 - 2.2. (56) 

The absence of long-range order has been traced back in these cases to domain wall pinning 
in metastable states 000. 

A short discussion of the domain relaxation ist found inta, which gives the following 
prediction for the non-equilibrium correlation length El 

R c (t)^^(c+^hx(t/t )) (57) 



h 2 \ r 

where C is a constant of order unity. This logarithmic growth has been seen experimen- 
tallyfi In order to explain the experimental data of diluted antiferromagnets in an exter- 
nal field at low temperatures, it is necessary to take into account also pinning by random 
bonds @. 



5 Miscellaneous 

With the present review I have tried to summarize some of the more recent activities in the 
theoretical investigation of the Ising model in a random field. Clearly, because of limitation 
of space and time not all new and interesting developments found their place in this review. 
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To mention one, Dahmen, Sethna and co-workers have studied the zero tempera- 
ture random field Ising model as a model for noise and avalanches in hysteretic systems. 
Changing the external field H(t) adiabatically from its initial state where all spins point 
downwards, at small disorder (i.e. for small h) the first spin to flip easily pushes over its 
neighbors, and the transition happens in one burst (an infinite avalanche). On the other 
hand, at large enough disorder the coupling between spins becomes negligible, and most 
spins flip by themselves, no infinite avalanche occurs. Tuning the amount of disorder in the 
system, i.e. the random field strength h, one finds a non-equilibrium critical point h c , H c 
where infinite avalanches disappear. At this point there is a universal scaling law for the 
magnetization m = M — M c {h c ) 

mn\h- h c fM((H - H c )/(h - h c f 5 ). (58) 

The critical exponents are believed to be those obtained for the equilibrium random 
field Ising model by dimensional reduction El. The relation between the critical field H c of 
Section 4.1 and H c is not clear at present. 

Other work not mentioned in this review considers the Ising model in a random field 
in d = 1 dimension or on special lattices, on which exact solutions are possible. I also left 
out the random field Ising model in a transverse field, the m — > oo limit of the random field 
Ising model etc.. We refer here the reader to the earlier review with P. Rujani or to the 
"Current Contents". 

The author of this review apologizes to all colleagues whose contribution has been 
omitted. He is conscious of its deficiencies, but nevertheless he hopes that it may be of 
some value for a first rough orientation in this field. 
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